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Abstract 

Here we present an algorithm to find elementary first integrals of rational second 
order ordinary differential equations (SOODEs). In [17], we have presented the first 
algorithmic way to deal with SOODEs, introducing the basis for the present work. 
In [18], the authors used these results and developed a method to deal with SOODEs 
and a classification of those. Our present algorithm is based on a much more solid 
theoretical basis (many theorems are presented) and covers a much broader family 
of SOODEs than before since we do not work with restricted ansatz. Furthermore, 
our present approach allows for an easy integrability analysis of SOODEs and much 
faster actual calculations. 



Keyword: Elementary first integrals, Second Order Differential Ordinary Equa- 
tions 

PACS: 02.30.Hq 

1 E-mails: javellar@dft.if.uerj.br, lduarte@dft.if.uerj.br, sduarte@dft.if.uerj.br and 

damota@dft.if.uerj.br 



1 Introduction 



The differential equations (DEs) are the most widespread way to formulate the 
evolution of any given system in many scientific areas. Therefore, for the last three 
centuries, much effort has been made in trying to solve them. 

Broadly speaking, we may divide the approaches to solving ODEs in the ones 
that classify the ODE and the ones that do not (classificatory and non-classificatory 
methods) . Up to the end of the nineteenth century, we only had many (unconnected) 
classificatory methods to try to deal with the solving of ODEs. Sophus Lie then 
introduced his method [1, 2, 3] that was meant to be general and try to solve any 
ODE, i.e., non-classificatory. Despite this appeal, the Lie approach had a short- 
coming: namely, in order to deal with the ODE, one has to know the symmetries 
of the given ODE. Unfortunately, this part of the procedure was not algorithmic 
(mind you that the classificatory approach is algorithmic by nature). So, for many 
decades, the Lie method was not put to much "practical" use since to "guess" the 
symmetries was considered to be as hard as guessing the solution to the ODE itself. 
In [4, 5], an attempt was made to make this searching for the symmetries to the 
ODE practical and, consequently, make the Lie method more used. 

So, despite all these efforts, a non-classificatory algorithmic approach was still 
missing. The first (semi) algorithmic approach applicable to solving first order 
ordinary differential equations (FOODEs) was made by M. Prelle and M. Singer [6]. 
The attractiveness of the PS method lies not only in the fact that it is based on a 
totally different theoretical point of view but, also in that, if the given FOODE has a 
solution in terms of elementary functions, the method guarantees that this solution 
will be found (though, in principle it can admittedly take an infinite amount of time 
to do so). The original PS method was built around a system of two autonomous 
FOODEs of the form x = P(x,y), y = V(x,y) with P and V polynomials in 
C[x,y] or, equivalently, the form y' = R(x,y), with R(x,y) a rational function of 
its arguments. 

The PS approach has its limitations, for instance, it deals only with rational 
FOODEs. But, since it is so powerful in many respects, it has generated many 
extensions [7, 8, 9, 10, 11, 12, 13, 14] 

Nevertheless, all these extensions deal only with FOODEs. In particular, the 
second order ordinary differential equations (SOODEs) play a very important role, 
for instance, in the physical sciences. So, with this in mind, we have produced [17] 
a PS-type approach to deal with SOODEs. This approach dealt with SOODEs that 
presented elementary 2 solutions (with two elementary first integrals). 

Here, we present an algorithm that comes supported by several theoretical re- 
sults. It proves to be very fast (in real applications) and, although being not com- 
pletely general, we could not find any interesting example where it could not be 
used. 

In section 2, we present the state of the art up to the present paper. In the 
following section, we introduce the algorithm to find the first integral. In section 4, 
we present some examples of SOODEs where the algorithm is successful. Finally, 
we present our conclusions and point out some directions to further our work. 

2 For a formal definition of elementary function, see [15]. 



2 



2 Earlier Results 

In the paper [6], one can find an important result that, translated to the case of 
SOODEs of the form 

// M(x,y,y') , 

y = "T77 k = nx, y, y ), (i) 

N(x,y,y') 

where M and N are polynomials in (x,y,y') 3 , can be stated as: 

Theorem 1: If the SOODE (1) has a first integral that can be written in terms 
of elementary functions, then it has one of the form: 

m 

I = w + ln(u)j), (2) 

i 

where m is an integer and the w's are algebraic functions 4 of (x,y,y'). 
The integrating factor for a SOODE of the form (1) is defined by: 

R^-y'') = d -^l (3) 
dx 

where represents the total derivative with respect to x. 

Bellow we will present some results and definitions (previously presented on [17]) 
that we will need. First let us remember that, on the solutions, dl = I x dx + I y dy + 
I y i dy' = 0. So, from equation (3), we have: 

R((f> dx - dy') = I x dx + I y dy + I y , dy = dl = 0. (4) 
Since y' dx = dy, we have 

R[((p + Sy')dx-Sdy-dy']=dI = 0, (5) 

adding the null term S y' dx — S dy, where S is a function of (x,y,y'). From equation 
(5), we have: 

I x = R(4> + Sy'), 

Iy = ~RS, (6) 
Iy' = —R, 

that must satisfy the compatibility conditions. Thus, defining the differential oper- 
ator D: 

D = d x + y'd y + ^d yl , (7) 
after a little algebra, that can be shown to be equivalent to: 



D[S] = -cj) y + S(P y , + S 2 , (8) 
D[R] = -R{S + <f> y ,), (9) 

Ry = Ry'S + Sy'R. (10) 



3 From now on, f denotes df/dx. 

4 For a formal definition of algebraic function, see [15]. 
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Combining (8) and (9) we obtain 

D[RS] = -R4> y . (11) 

We will now briefly explain the algorithm proposed on [17]. In that paper, we 
made the conjecture that it is always possible to have an first integral I such that 
RS is a rational function of (x,y,y'), if the corresponding SOODE presented an 
elementary solution. From equation (11), we then see that this conjecture implies 
that both R and S are also rational functions. 

Using this into equation (8) and trying to solve it, we are left with the task of 
solving a third degree algebraic systems of equations on the desired coefficients of 
the polynomials defining the rational function S. Supposing that this stage was 
overcome, we would then substitute the S just found in equation (9) and apply the 
a Prelle-Singer type approach to find R. 

Basically, that was the novelty presented on [17], as far as we know, the first 
semi-algorithmic approach to tackle SOODEs. 



3 The Algorithm 

In the present section, we will present a non-classificatory algorithm to search el- 
ementary first integrals for SOODEs. The method is based on a Darboux type 
procedure and, analogously to the Prelle-Singer [6] approach for FOODEs, it guar- 
antees that, for a class of SOODEs, it will eventually find the desired first integrals. 
The algorithm relies on several theoretical results that will be introduced on the 
next sub-section. 



3.1 Theoretical Foundations 

Let us start then by a corollary to theorem 1 concerning S and R. 

Corollary 1: If a SOODE of the form (1) has a first order elementary first 
integral then the integrating factor R for such a SOODE and the function S defined 
in the previous section can be written as algebraic functions of (x,y,y'). 



Proof: Using the above mentioned result by Prelle and Singer, there is always 
a first integral I = wq + J2T c i In^i) f° r the SOODE. So we have, using equation 

(3), 

R{ Jl ~ V " ] = h + V ' Iy + V " Iy ' ^ R = ~ Iy ' (12) 
where I u = d u I. From equation (2), we have: 



l y , = w 0y , +Y,a-^j L - (13) 

<JJ 



Then I y i is an algebraic function of (x,y,y') and, by equation (12), so is R. 
From equations (6), one can see that: 

-r i \-^7Tl w iy 

S=±L= W ° y + ^ d4) 

T , i \-^m w iy' v ' 
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Therefore, S is also an algebraic function of (x,y,y').0 

In order to produce the further mathematical results we need, we are going to 
re-define the functions R and S. 

Suppose we have a SOODE of the form (1). We will then define the integrating 
factor as: 

K(M - Ny") = d/(X ^' y/) = Ix dx + I y dy + /„, dy' = 0. (15) 
This is equivalent to use the following re-definition: 

n=§. (16) 

Proceeding analogously to what we did with equation (4), we get: 

IZ[(M + Sy')dx-Sdy-Ndy']=dl = 0, (17) 

adding the null term Sy' dx — S dy, where S is a function of (x, y, y'). 
This is equivalent to use the following re-definition: 

S = SN. (18) 

Using equations (16) and (18) into equations (9), (10) and (11), we get: 

vim 

_J_i = _ {S + N x + N y y' + M y! ), (19) 

lZ y N - TZNy = N 2 (Jl yl S + S y > K) (20) 

S^ + V[S\ = —(N My — M Ny) (21) 

where V is defined as V = N D. Since TV is polynomial and R and S are algebraic, 
so are 1Z and S. 

Now, we will introduce some theorems that will be the basis of our algorithm. 

Theorem 2: Consider a SOODE of the form (1), that presents an elementary 
first integral I. If S is a rational function of (x, y, y'), then the integrating factor TZ 
for this SOODE can be written as: 



K = ftp? (22) 



where pi are irreducible polynomials in (x, y, y') and n{ are non-zero rational num- 
bers. 



To prove Theorem 2 we will need the following lemma: 

Lemma 1: Consider a function F of (x,y,z). If the differential of F can 
be written as dF = A(X dx + Y dy + Z dz), where X, Y and Z are polynomial 
functions of (x, y, z) and A is an algebraic function of (x, y, z), then the first order 
ordinary differential equation defined as 

^ - -- (23) 
dx~ Y {16) 
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where z is regarded as a parameter, has F(x,y,z) = C (where C is a constant) as 
a general solution 

Proof of Lemma 1: Consider the first order ordinary differential equation denned 
by (23). If / is a function of (x,y) such that 

(Yd x -Xd y )[f(x,y)}=0, (24) 

then f(x,y) = k (where k is a constant) is a general solution of (23). Applying 
(Y d x — X d y ) to F(x, y, z) we get Y F x — X F y . But, by hypothesis, F x = A X and 
F y = AY leading to 

(Yd x - Xd y ) [F{x,y,z)\ = Y AX - X AY = 0. 

This implies that F(x,y,z) = C is a general solution of (23). □ 

Proof of Theorem 2: If the hypothesis of the theorem are satisfied, S is a rational 
function of (x,y,y'). So we can write S = P/Q where P and Q are polynomials. 
Substituting this into equation (17) we get 

^ [{MQ + Py')dx-Pdy- NQdy'] = dl = 0, (25) 

By using Lemma 1 we have that I(x, y, y') = C is a general solution of the FOODE 
defined by 

dy = MQ + Py' 

dx P { ' 

If the polynomials (M Q + P y') and P have a common factor we can write (M Q + 
Py') = Ti tu , P = T 2 tn and (26) as 

T- = ( 27 ) 
dx T 2 y ' 

Since I(x, y, y') is an elementary function (by hypothesis) then, by the theorem of 
Prelle and Singer [6], there exists an integrating factor R\2 for the FOODE (27) of 
the form 

Rn = iifr (28) 

i 

where /j are irreducible polynomials in (x, y) and m; are non-zero rational numbers. 
By other side, from equation (25) we have that 

Ix = -^t 12 T 1 ; I y = ^t 12 T 2 , (29) 
and, so, ^ t± 2 is an integrating factor for (27). This implies that 

Rl2 = HI) ^12, (30) 

where is a function of the first integral /. From (30) is easily seen that 

F(I)n = R 12 Q- (31) 

tl2 
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and we may notice (see (28)) that J 7 (I) H can be expressed as n p"} % ■ But, since 1Z is 
an integrating factor for the SOODE (1), so is R = ^(1)71. Finally, the reasoning 
we apply to the pair of variables (x,y) when we make use of Lemma 1, could be 
carried out analogously with the pairs (x, y') and (y, y') as well. So, we can conclude 
that there exists an integrating factor of the form (22). □ 



Theorem 3: Consider a SOODE of the form (1), that presents an elemen- 
tary first integral I. If S (defined above) is a rational function of (x,y,y f ) (S = 
P/Q, where P and Q are polynomials in (x,y,y')), then (1Z/Q) | V[TZ/Q], i.e. 
V[1Z/Q]/(1Z/Q) is a polynomial. 

Proof of Theorem 3: Substituting S = P/Q in the compatibility conditions 
(19,21), we get: 

Villi 

Q^ = -P-Q(N X + N y y' + M yl ), (32) 

r_m + 9m^pm =_ (NMy _ MNt) (3S) 

Multiplying (33) by Q, one gets: 

P ^ + V[P] - P ^ = -Q (N M y - M N y ) (34) 

and finally 

p ^ _ = _ v[p] _ Q[NMy _ MNyY (35) 

If one adds — V[Q] to both sides of (32), one gets: 

Q - ^) = -V[Q\ -P-Q(N X + Nyy' + M y ,) (36) 

Re- writing (36,35), we have: 

Q (^r) = ~ V [Q\ -P-Q{N X + N y y' + M y/ ) (37) 

and 

" D ^) = -V[P] -Q(NMy-MNy). (38) 

Since, by definition, P, Q, M and N are polynomial and V is a differential opera- 
tor with polynomial coefficients, the right-hand side of both (37,38) are polynomial. 
Therefore, in principle, ^J^q^ is rational. So, let us represent it as A/B, where A 
and B are polynomial that do not have any common factors. By doing that, one 
can write (37,38) as: 

Q (|) = Vi (39) 

A 
B 
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where Vi = -V[Q] - P - Q (N x + N y y' + M y ,) and V 2 = —V[P] — Q (N M y — M N y ). 

In order to satisfy (39,40) simultaneously, since A and B do not have common 
factors, it would be necessary that B\Q and B\P. But, again by hypothesis, P 
and Q do not have common factors. So, one can conclude that we have B = 1. 
Therefore, A/B = V ^/q ] = polynomial. □ 

Corollary 2: Consider a SOODE of the form (1), that presents an elementary 
first integral I. If S is a rational function of (x,y,y') (S = P/Q, where P and Q 
are polynomials in (x,y,y')), then 1Z/Q can be written as: 

K/Q = l[vr (41) 

i 

where Vi are irreducible eigenpolynomials (in (x,y,y')) of the V operator and m,i 
are non-zero rational numbers. 



Proof of Corollary 2: From theorem 2 we have that 1Z = WiP™ 1 , where pi are 
irreducible polynomials in (x, y, y') and are non-zero rational numbers. Since Q 

k 

is polynomial, we have Q = Ylj , where qj are irreducible polynomials in (x, y,y) 
and kj are non-zero positive integers. So, 

n /Q = Y[v^ (42) 

i 

where Vi are irreducible polynomials (in (x,y,y')) and m; are non-zero rational 
numbers. Since, by theorem 3, V[R/Q]/(R/Q) is polynomial, we have: 

So, since the Vj's are irreducible, t>i|D[t>;].n 



3.2 The Algorithm Itself 

In this section, we will make use of the mathematics constructed above in order to 
produce a "semi- algorithm" to deal with SOODEs of the class defined in (1). 
From theorem 3 we can write: 

n/Q = T = => ^ = Em, ^1 = Dmift (44) 

where Vi are irreducible eigenpolynomials (in (x,y,y')) of the D operator and the 
g' Si are the corresponding cofactors. 

From equations (37,38), we can then write: 

Q (^P) = Q^2m^ = -V[Q]-P-Q(N x + N y y' + M y/ ) (45) 

and 

P (^-^=PY, m *9 i = -V[P]-Q(NM y -MN y ). (46) 

Equations (45,46) will be the basis of our procedure. Let us begin discussing our 
procedure by talking us through the main steps before formalizing the algorithm. 
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As it will become clear as we go along, the first step of our procedure is the 
most costly one, i.e., time consuming. In the same way as in the Prelle-Singer 
procedure for FOODEs [6], determining the eigenpolynomials and corresponding 
cofactors for the T> operator is a task that grows as the degree of such polynomials 
grow. Our present procedure start by performing this search. Assuming that this 
search succeeded and we found V{ and gi for some degree deg, checking (45,46), 
we see that it remains to be determined the following set of variables: {mj, P, Q}, 
where nii is a constant and P and Q are polynomials. How to determine these? 
Again, by inspecting (45,46), we note that if we construct a generic polynomial 
Q (in {x,y,y'))) of some degree degQ, we can infer the maximum degree of the 
polynomial P (degp). So, we construct such polynomials (the most general ones for 
the corresponding degrees) and solve (45,46) for the coefficients of such polynomials 
and for rrij. By doing that, we would have found mj, P and Q and, consequently, 
T = YliV™*. Since we have 1Z/Q = T, we would have found R, the integrating 
factor for the SOODE in question! 

Once R and 5 5 have been determined, we have all the partial first derivatives 
of the first order differential invariant, I(x, y, y'), which is constant on the solutions 
(see equations (6)). This invariant can then be obtained as 



Basically, there are a few scenarios that can arise in the case of SOODEs (in regard 
to the first order differential invariants): One can find two pair of independent S 
and pleading to two independent I's. In this case, we would have fully integrated 
the SOODE (by solving for y', using one of the invariants, and substituting in the 
other to find y). On the other hand, it is possible that we can only find one such 
invariant and, in that case, we can only reduce the SOODE to a FOODE. 

This overview above indicates the main trust of our procedure, bellow we will 
present a step-by-step modus operandi of how to implement the above scenario. 

• Steps of the Algorithm 

1. Set Deg = 1. 

2. For the given SOODE, we determine the eigenpolynomials and the asso- 
ciated cofactors, up to degree Deg. 

3. Set DegQ = 1. 

4. Construct a generic polynomial Q (in {x,y,y')) of degree DegQ and a 
generic polynomial P (in (x, y, y')) of degree Degp = DegQ+M.AX.(degM — 



5. Try to solve equations (45,46) for the coefficients defining Q and P and 
for 771 j . 

6. If we are successful, go to step 7. In the opposite case, if DegQ < 10 Deg, 
we make DeqQ = DegQ + 1 and return to step 4. If DegQ = 10 Deg, set 
Deg = Deg + 1 and return to step 2. 




l{ R+ w[S m+Sv ' )ix -S 




l,deg N ) 6 . 



5 Please remember that S = P/Q. 

6 Where degu and degw are the degree of M and N, n 



•espectively. 
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7. Now that we have S and 1Z we, after checking if that solution satisfy 
equation (20), in the affirmative case, using equation (47), calculate the 
associated first order differential invariant. If the solution does not satisfy 
(20), set DeqQ = DegQ + 1 and return to step 4 above. 

4 Examples 

In the present section, we are going to show some examples of the usage of our 
method. In the first example, we are going to point out the improvements that have 
been made in comparison to the first algorithmic approach to solving SOODEs pre- 
sented on [17]. Later we will apply the method to another two potentially interesting 
physical examples found on a very interesting paper [18] where the authors use our 
results [17] and construct a classification of SOODEs and present different heuristics 
to integrate each one of them. We will comment on that classification in the light 
of our present theoretical results and issuing algorithm and make some remarks 
about the capability of our method to analyze the integrability of the corresponding 
SOODEs. Finally, on a more mathematical tone, we will show that the method can 
solve equations previously unsolved by other powerful techniques. 

4.1 First Example 

A rich source of non-linear DEs in physics are the highly non-linear equations of 
General Relativity. Einstein's equations are, of course, in general, partial DEs, 
but there exist classes of space-times where the symmetry imposed reduces these 
equations to ODEs in one independent variable. One such class is that of static, 
spherically symmetric space-times, which depend only on the radial variable, r. 
The metric for a general statically spherically space-time has two free functions, 
A(r) and fi(r) say. On imposing the condition that the matter in the spacetime is 
a perfect fluid, Einstein's equations reduce to two coupled ODEs for A(r) and n(r). 
Specifying one of these functions reduces the problem to solving an ODE (of first 
or second order) for the other. 

Following this procedure, Buchdahl [19] obtained an exact solution for a rela- 
tivistic fluid sphere by considering the so-called isotropic metric 

ds 2 = (1 - /) 2 (1 + /)- 2 dt 2 - (1 + /) 4 [dr 2 + r 2 (d# 2 + sin 2 6 d<p 2 )} 

with f = f(r). The field equations for f(r) reduce to 

//" " 3/' 2 - r~ l ff = 0. 

Changing to the notation of this paper, with y(x) = f(r), we get: 

yy"-3y' 2 -(l/x)yy' = 0. (48) 

Let us now apply our new method (section 3) to the above equation. First, we 
have to calculate the eigenpolynomials of degree=l (and the associated cofactors) 
for the differential operator associated with equation (48) given by: 

V = x y 8 X + x y y' 8 y + y' (3 y'x + y) d y , . (49) 
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These are found to be: 



vi = x 9i = y (50) 

V2 = V 92 = xy' (51) 

vs = y' 53 = 3 y' x + y (52) 

Now, following the steps of the algorithm (given on section 3), we can find two 
independent solutions: 



• First solution 



S = -3y'x (53) 

R = ih (54) 



Using equation (47), we find, for this pair of S and TZ: 



• Second solution 



i = 4r- (55) 

y 6 x 



_ x(2yx + 3y' + 3x 2 y') 

s ~ TT^ (56) 

Using equation (47), we find, for this pair of S and 1Z: 

I = y' + *y + y\ (58) 

y 6 x 

In comparison to the method we have presented on [17], we see that our present 
approach has some remarkable differences: It covers a much broader "universe" of 
SOODEs since, in the present algorithm, we do not have restrictions imposed on 
1Z (in the algorithm presented on [17], we made a conjecture that led to 1Z being 
necessarily rational). The present approach is much more theoretically sound than 
the previous one, we produced a lot of theorems to support the algorithm. From 
a practical point of view, this new approach is much more efficient. For example, 
equation (8) will generated a system of third degree algebraic equations in the 
coefficients we want to determine. On the other hand, in the algorithm introduced 
in section 3, equations (45, 46) will mostly generate first degree algebraic equations 
since only the terms (Q J2 m i9i) an d (P J2 m i9i) produce second degree equations 
on the desired coefficients. 

In [18], based on our results [17], the authors produced an way of finding a 
second independent solution for the S and 1Z from the first one 7 . As it is clear 
from the above, our present approach finds both pairs independently. They also 
classify the SOODEs they dealt with into three types. The example we presented 
above (that was originally treated in [17]) was also used by the authors in [18] as a 
representative of type I SOODEs. 

In the following two sub-sections, we will introduce an example for each of the 
other two types. 



7 At the time they presented their paper, they used our results from [17] where we introduced the 
theory and were concerned with the finding of only one pair of S and 1Z 
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4.2 Second Example 

Now, we are going to consider the Helmholtz oscillator with friction: 

y" + c iy ' + c 2 y- (3y 2 = 0. (59) 

It is an well known result that the above equation is integrable for the choice 
C2 = IE c ! 2 P^]- As we shall see, our method takes care of such concerns, i.e., since 
the algorithm translates the solving of the SOODE into solving systems of algebraic 
equations, if we consider the arbitrary constants appearing on the SOODE (ci, c 2 , (3) 
as variables (as well as the coefficients defining the eigenpolynomials and cofactors), 
the solution of such algebraic systems will provide the regions of integrability nat- 
urally. In order to clarify this, let us then follow our procedure in more detail than 
in the previous example: 

In order to follow our algorithm, the first thing to do is to determine the eigen- 
polynomials (and corresponding cofactors) for the V operator: 

V = d x + y'd y + (-ci y'-c 2 y + (3 y 2 ) d y , (60) 

The equation to be solved is: 

V[v]=gv (61) 

where v is the eigenpolynomial and g is the cofactor. 

It turns out that for degree 1 and 2 there is no solution to (61). A generic 
polynomial of degree 3 is given by: 

v = ai + a 2 y 2 y' + a 3 yy' 2 + a 4 x 3 + a 5 y 3 + a 6 y' 3 + a 7 xyy' + a 8 x + 
ag V + aio x 2 + an yy' + a 12 y 2 + a i3 xy' + a u xy + a i5 y' 2 + 

ai6 x 2 y' + an x 2 y + a 18 xy 2 + a 19 xy' 2 + a 2 o y (62) 

Analyzing the T> operator (eq. (60)) and equation (61), one may conclude that the 
maximum degree for g is 1. So, 

g = h + b 2 x + 6 3 y + 64 y'. (63) 

Solving equation (61) will be translated into solving the following algebraic system: 



sys = (f3a 2 - a 5 63 = 0, [5 ai 3 - c 2 a 7 - a\ 2 b 2 - 


ais 01 - a i4 63 


= 0, 


-a 5 b 2 - a 18 63 + (3 a 7 = 0, 2 a w - a 7 b 3 - 


- a 2 b 2 - ais h 


= 0, 


-aio o 4 - ai 6 61 - Ci ai 6 - ai 3 b 2 + an = 0, - 


019 b 2 — a±Q 64 


= 0, 


-ai 6 63 - 


- an &4 — 07 6 2 


= 0, 


-aio 63 - ai7 h - au b 2 - c 2 a w = 0, -a w b x 


- a 8 6 2 + 3 a 4 


= 0, 


-a 4 6 2 = 0, -an b 2 - 04 63 = 0, -ai b 2 - 


- a 8 h + 2 aio 


= 0, 


-aio &2 - a 4 61 = 0, -ai 6 6 2 - a 4 64 = 0, 


-a 6 63 - a 3 64 


= 0, 


-ais &2 + /? ai 6 - ai 7 b 3 = 0, -a 3 6 2 - 


- a 7 64 - a i9 63 


= 0, 


-a 3 63 + 3 a 6 - a 2 64 = 0, -a 9 63 - c 2 an + ai 8 - 


a i2 61 + /?a 20 


= 0, 


-c 2 a 2 + (3 an - 


- ai2 o 3 - a 5 61 


= 0, 


-c 2 a 20 + a i4 - ai 63 - a 9 61 = 0, -3 cj a 6 + a 3 - 


• a i5 64 - a 6 61 


= 0, 
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2 P a 3 - a 2 63 - a 5 64 = 0, aig - 2 cj a i5 + an - a 20 64 - »15 &i = 0, 
— ae 64 = 0, 013 — ai 64 — 020 »i — a 2 o + a 9 = 0, 
a 8 - ai 61 = 0, -3 c 2 a 6 + 2 a 2 - 2 a 3 - a i5 fe 3 - a n 64 - a 3 61 = 0, 
-2 c 2 ai 5 + 2 012 - a 20 o 3 + a 7 - a 9 64 - an 61 - c; an = 0, 
-2 aig - a i3 64 - aig 61 - a 15 b 2 + a 7 = 0, 
-a 20 6 2 - a 8 64 - ai 3 61 + 2 a i6 - cj a i3 + a M = 0, -a 6 6 2 - aig 64 = 0, 
-2 c 2 aig - an b 2 - cj a 7 — a i3 6 3 + 2 ai 8 - a i4 64 - a 7 61 = 0, 
2 ai7 - c 2 ai 3 - a 9 6 2 - a i4 61 - a 8 6 3 = 0, 
-ai2 o 4 - a 2 61 - an 6 3 - cj a 2 - 2 eg a 3 + 3 a 5 + 2 /? ai 5 = 0) (64) 

A solution to the system above is: 

a 20 = 0,a 4 = -^/3an,a 9 = Ci 2 an, c 2 = ^ c 2 ,a^ = 0, 

ai 8 = 0, a i2 = 0, a 5 = 0, a i5 = 0, ai 6 = 0, 6 3 = 0, 64 = 0, a i9 = 0, 

6 

a 7 = 0, a 2 = 0, a 8 = 0, b 2 = 0, ai = 0, b\ =--ci, an = an, a i3 = 0, 

5 

4 

aio = - ci an, 017 = 0, 014 = 0, a§ = 0. (65) 



Leading to: 



V = 


-1"* 


9 = 


6 




6 2 


c 2 = 


25 Cj 



5 



4 2 2 4 . /2 

y +-c iyy +y 



(66) 

As we have mentioned, the solution via our procedure takes care of the analysis 
of integrability of the SOODE. The above value for c 2 is an well known result for 
the integrability of the Helmholtz system. 

So, following the remaining steps of our algorithm, solving equations (45) and 
(46) to find P and Q, we can get two independent solutions: 

• First solution 

s -25py 2 +Ac 1 2 y + lQc 1 y' 

10c J y + 25y / " 1 ) 

n _ 2 Cl y + 5 y' 

-50 y 3 + 12 Cl 2 y 2 + 60 a y y' + 75 y' 2 

Using equation (47), we find, for this pair of S and 1Z: 

I = 6 Cl x + 5 In (-12 Cl V + 50 p y 3 - 60 cj y y' - 75 y' 2 ) (69) 

• Second solution 

_ 6 Cl 2 y + 25 cj y'-25Py 2 

S ~ 25? (70) 

/ 9 4 4 \ — 5/6 

ft = y[-^(3y 3 + ^c 1 2 y 2 + -c 1 yy' + y' 2 ) (71) 
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Using equation (47), we find, for this pair of S and 7Z: 

i = I — (72) 

What new properties of our algorithm can be emphasized by this example? 
As mentioned, the method is blind as for the origin of the parameters we want to 
determine. So, one can use this to include the integrability analysis in a very natural 
way. We see that, since our method is based on very sound theoretical results (see 
section 3) about the nature of 1Z, we are left with a truly algorithmic procedure of 
the PS- type, i.e., we reduce the solving of a SOODE to an algebraic problem. 



4.3 Third Example 

The next example is also physically motivated, the force free Duffing- van der Pol 
oscillator: 

y" +(a + f3y 2 ) y' - iy + y 3 = Q (73) 

For this case, there are also some integrability aspects. We are not going to 
dwell on that now. This was already done on the last example. It is well known 
(see for instance, [21]) that the integrability occurs for 7 = —3 (3~ 2 , a = 4 f3~ l . 

Applying our method (section 3) to that equation, we first have to calculate the 
eigenpolynomials to the appropriate degree (and the associated cofactors) for the 
differential operator associated with equation (73) given by: 



These are: 



V = 8 X + y'd y + (-y 1 a- Py 2 y' + iy -y 3 ) d y , (74) 

n = — ^- = — (75) 

(3j/ + /?V + 3j//3) 3 
V2 = 92 = — ^ (76) 

Now, following the steps of the algorithm (given on section 3), we can find two 
independent solutions: 



• First solution 



2 y 2 + 1 

S = (77) 
11 = (3y + /?V + 3 y'/?) (7§) 



Using equation (47), we find, for this pair of S and 7Z: 

3x + ln(3y(3a-9y + [3 2 y 3 + 3y'[3) [3 



(79) 



Second solution 



s = 4y'(3 + y'^y 2 + 3y + y 3 (3 2 (gQ) 
K = V == (81) 
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Using equation (47), we find, for this pair of S and 1Z: 
I = -(Sy + y 3 p 2 + 3 y'p) ^ ^ + 2 In ( foy + y*p + Zy'f3 _ y ^2/3^ _ 



In ( (3 y + y 3 (3 2 + 3 y'(3) ^ + ^3y + y*P + 3y'/? j, /3 2 / 3 + y 2 /3 4/3 ) + 



2V3arctanf^ 3^/3 "J ^ 

In this example, we skipped the integrability analysis (that was already exem- 
plified on the previous one but could be done in the same fashion here) and followed 
the same procedure as in the two first examples to find the functions S and 1Z. 



4.4 Some Comments on the Examples Above 

Before we present the fourth (and final) example, let us make some considerations 
on the three previous ones: In [17], we have introduced the basis for an approach 
to algorithmically solve SOODEs. Later on, [18] picked up some of these ideas and 
have developed a different approach to the task. They also used our function S as 
one of the main ingredients of their method. Shortly, they say they use the same 
ansatz as we did in [17] for S: 

s = a(x,y) + b(x,y)y' 
c(x,y) + d(x,y) y' 

Actually, this is a more general form than the one we considered in [17]. As we 
do here, we consider S to be a plain rational function on (x, y, y') (not rational in 
respect only to y'). Using their format (83), our case would correspond to (a, b, c, d) 
being polynomials on (x,y) 8 . This is an important point in understanding the 
differences between our present algorithmic method and the approach presented on 
[18]. To further our discussion, let us introduce the results for S ad R, presented 
on [18], for the three examples above: 

• For example 1 above (4.1) 

The authors of [18] have found two pairs of S and R: 

Si = — R\ - 



x y 3 x 



v' 1 

S 2 = ~- R2 = ^~ (84) 
x y° x 



For example 2 above (4.2) 

For this case, the authors have found: 

,2 1 4c ? „ , 2ci „,/ 



y + y' 



s 2 = ^!±^l±li R2 = _y> e ci - (85) 



v 



8 Please note that, in our considerations, we were not limited to numerators and denominators of 
degree 1 in y', etc. 



15 



• For example 3 above (4.3) 

In this third example, the authors of [18] have found only one pair of S and 
R: 

Si = ^ + Py 2 R! = eT- (86) 

In order to obtain the above results, the authors also used an ansatz for R: 

R = A{x 1 y) + B{x 1 y)y l . (87) 

They use the two ansatz (for S and R) into equations (8), (9) and (10) thus ob- 
taining a system of partial differential equations for the unknown functions (a,b,c,d, 
A,B) of (x,y), this can be very trick and not algorithmic in essence (see example 
3 above where the authors could only find one independent pair of S and R). On 
the other hand, our method only consider that S has to be a rational function of 
(x, y, y') and, due to the theoretical results presented on section 3, we know the gen- 
eral form of R and that led to an algorithmic Darboux type procedure to determine 
S and R. 

As can be seen from equations (85, 86), the method presented on [18] can find, in 
principle, non-algebraic S and R (our algorithm deals with rational S and algebraic 
R). But, for the examples presented on [18], the authors found only rational S 
functions and, therefore, in principle, inside the applicability of our method. As 
mentioned already, the point is that our approach is algorithmic (in the Prelle-Singer 
sense) and, computationally speaking, much faster to apply. Basically, we convert 
solving a SOODE to solving an algebraic system while the method presented on 
[18] transforms the solving of the SOODE in solving a system of partial differential 
equations. 

The three examples above are physically motivated examples and were previously 
dealt with in [17] (example 1) and [18] (examples 1,2,3). In this latter reference, 
these were classified into three different types. It is worth to point out that, in 
our algorithm, that is not necessary, i.e., all examples are treated equally since our 
procedure is non-classificatory. 

In order to conclude this section with examples, we will pass now to a more 
mathematical note and present an academic SOODE with the intention to show 
that the present approach is capable of solving it while many powerful techniques 
fail. 



4.5 Fourth Example 

Let us consider the following SOODE: 

y" = fr^ ^ 

-y - 1 + 3xy' 

Let us now apply our new method (section 3) to the above equation. First, we 
have to calculate the eigenpolynomials of degree=l (and the associated cofactors) 
for the differential operator associated with equation (48) given by: 

T> = (-y - 1 + Sxy') d x + {-y'y - y' + 3y' 2 xj d y - y' 2 d y >. 
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These are found to be (up to degree=l): 

vi =y' gi = -y' (90) 

Now, following the steps of the algorithm (given on section 3), we can find two 
independent solutions: 

• First solution 

_ y'{-2y + 3xy>-2) 

(- y -l + 3xy') ^ 
K = -y-l + 3xy' (92) 

Using equation (47), we find, for this pair of S and 1Z: 

I = 2x y' 2 y + 2 y' 2 x - 2 y /3 x 2 - | yV -^yy'-^y'. (93) 

• Second solution 

U = y' [ * + ~> (95) 
Using equation (47), we find, for this pair of S and 1Z: 



I = arcsin 



yVs + Vs 



{J(6xy'-3y-3) 2 + {yV3 + V3) J 



+ ^ In (y>) . (96) 



The great advantage of having a non-classificatory, algorithmic approach is that 
we can tackle every SOODE in the same manner. The above example, despite 
its simple appearance, could not even be reduced by the powerful Maple solver 
(release 9.5). In our method, the most costly part of the algorithm is to determine 
the eigenpolynomials (and corresponding cofactors) and, for this example, this is 
actually very simple. 

5 Conclusion 

In [17], we have developed a method, based on a conjecture, to deal with SOODEs 
that presented an elementary solution (possessing two elementary first integrals). 

In that same paper, we have introduced a function S to transform the Pfaffian 
equation related to the particular SOODE under consideration into a 1-form pro- 
portional to the differential of the first integral. That function S was instrumental 
in finding the integrating factor for the SOODE. 

Here, in the present paper, we introduce many theoretical results concerning 
that function S and present an way to calculated it and the integrating factor R via 
a Darboux-type procedure. 

Briefly, we construct a differential operator V, 

V = N D = Nd x + y' Nd y + Md y ,, (97) 
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extracted from the SOODE, and the corresponding eigenpolynomials and cofactors 
will be the building blocks of S and R. Furthermore, the procedure is semi- algo- 
rithmic and, given enough time, if the solution exists, it will find it (as does the 
Prelle-Singer approach for FOODEs). 

From an operational point of view, our method is very sound, since, as mentioned 
in sub-section 4.2, our algorithm converts the solving of the SOODE into solving 
(essentially) first degree algebraic equations. 

As a consequence of this, our approach is capable of analyzing the integrability 
regions for the SOODE (for the case where it presents undetermined parameters). 
Since we translate the SOODE into an algebraic system, if we consider those pa- 
rameters as variables, the method solves the system given the values for which there 
is an integration possible. 

Our method is not general since it is limited to the cases where there exists 
an elementary first integral and, furthermore, to the cases where S is rational. 
Actually, in respect to the first limitation mentioned above, it is possible that our 
method solves a case where the first order differential invariant is not elementary 
(see example 2, sub-section 4.2). The point is that, if S and R are of the right format 
(i.e., rational S and R = JJiP^, where pi are irreducible polynomials in (x,y,y') 
and rii are non-zero rational numbers.), our method will work. The method is sure 
to work if there is an elementary first order invariant. We would like to elaborate a 
little bit bout the second restriction mentioned above (S being rational): although 
all of our theoretical results apply to that particular situation, we could not, so far, 
find any an interesting (not easily solvable) example where S is an algebraic function 
of (x, y, y') (the general case for elementary first order invariants, see Corollary 1 
above). So, the restriction does not seem to be very restrictive. For instance, in 
[18], all cases of physical interest are covered by our method, i.e., present rational 
S. 

In regard to future work, many extensions of the present paper can be pur- 
sued: We intend to further our work to include Liouvillian first order differential 
invariants (along the lines we followed, in dealing with FOODEs [13]). We can also 
include SOODEs with elementary functions (see [7, 14], for related work applied to 
FOODEs). As hinted above, we also intend to further analyze the capabilities of 
our method in the determination of integrability regions for the parameters present 
on the SOODE. 
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